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The ground states of few electrons confined in two vertically coupled quantum rings in the presence 
of an external magnetic field are studied systematically within the current spin-density functional 
theory. Electron-electron interactions combined with inter-ring tunneling affects the electronic struc- 
ture and the persistent current. For small values of the external magnetic field, we recover the zero 
magnetic field molecular quantum ring ground state configurations. Increasing the magnetic field 
many angular momentum, spin, and iso-spin transitions are predicted to occur in the ground state. 
We show that these transitions follow certain rules, which are governed by the parity of the number 
of electrons, the single particle picture, the Hund's rules and many-body effects. 
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I. INTRODUCTION 

The physics of semiconductor nanostructures has been 
a subject extensively studied since the experimental re- 
alization of quantum dots in the 1980s 1 -. This system is 
very interesting due to their similarity with atoms and 
the facility in controlling their electronic, magnetic and 
optical properties 2 -. With the advances in experimen- 
tal techniques to fabricate nanostructures, a novel ring- 
shaped nanostructure was pursued and realized through 
several different approaches, e.g. nano-lithography (e.g. 
atomic force microscope patterning) 3 -, droplet MBE 
epitaxy^ and strain induced self-organization^. These 
ring-shaped nanostructures, the so-called quantum rings 
(QRs), are known for their Aharonov-Bohm effect and 
its persistent current^ where the quantum interference 
phenomenon leads to oscillations in the current. 

Such ring-shaped nanostructures can also be cou- 
pled in the lateral or vertical configuration forming a 
"benzene-like" artificial molecule. When QRs struc- 
tures are coupled they open the avenue of controlling 
and manipulating some fundamental quantities, as e.g., 
the electron-electron interaction^ ^ 10 ! 11 ] 12 ! 13 , the in- 
teger and fractional Aharonov-Bohm oscillation^, the 
electron relaxation-^ and the coupling of direct-indirect 
excitonsi^. Quantum ring molecules (QRMs) were 
synthesized experimentally using MBE technology in 
the form of vertically stacked layers of self-assembled 
QR a 17 i 18 and concentric double QR a 19 i 20 . And recently, 
the Aharonov-Bohm oscillation was observed in self- 
assembled InAs/GaAs quantum rings containing only a 
single confined electron 2 ^. 

In a previous work 2 ^, we investigated theoretically 
the persistent current in two vertically coupled quantum 
rings (CQRs) containing 6 electrons in a wide range of 
inter-ring distances up to ~ 10 4 A. The motivation of this 



'Present address: Department of Physics, University of Cal- 
ifornia San Diego, La Jolla, California 92093-0319, USA, 
email:lcastelano@physics. ucsd.edu 



work was to understand the inter-ring quantum tunnel- 
ing effects on the persistent current of two interacting 
coupled rings. In order to analyze such effects, we con- 
sidered two different situations. First, we assumed that 
each quantum ring (QR) contains 3 electrons and inter- 
act with each other only via Coulomb potential. In the 
second situation, we included the possibility of tunneling 
between electrons localized in different QRs. We com- 
pared both situations and we found that the persistent 
current is altered significantly by the quantum tunneling, 
which allows the exchange interaction between electrons 
localized in different QRs. Moreover, we found that an 
applied vertical gate can be used to control the persistent 
current in such a system. 

In the present work, we apply the current spin-density 
functional theory (CSDFT) to determine systematically 
the ground states of two vertically CQRs in the presence 
of an external magnetic field. The two quantum rings are 
tunnel coupled and form a simple artificial molecule. Ear- 
lier, this method was employed to investigate the ground 
state properties of a single quantum rin g 23 ! 24 and it was 
proved to be a useful tool to determine the properties of 
such systems where confinement, Coulomb interaction, 
spin polarization, and magnetic field are present at the 
same time. In the two CQRs case, the inter-ring coupling 
plays an important role and by varying the inter-ring dis- 
tance between the QRs and the applied magnetic field 
we found a rich variety of ground state configurations for 
the systems of 7V=3,4,5, and 6 electrons. Such results 
are compiled in "phases diagrams" . Increasing the in- 
terring distance, the isospin quantum number decreases 
monotonously and the interring exchange-correlation in- 
teraction plays an essential role forming new molecular 
states. Also, we found some rules based on the single- 
particle picture to estimate the ground state configu- 
ration, which might be used as a guide to experimen- 
tal analysis in the future. The persistent current as a 
function of either magnetic field or distance between the 
QRs is determined for a different number of confined 
electrons iV=3,4,5, and 6, thereby completing our pre- 
vious work^ 2 .. We verify that the ground state configura- 
tion and the persistent current are inter-twined. There- 
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fore, through measurements of the persistent current in 
CQRs, the ground state configuration can be accessed 
experimentally. Furthermore, we evaluate the persistent 
spin-current, which is given by the difference between the 
spin-up current and spin-down current. We verify that 
the persistent spin-current is closely related to the ground 
state configuration of the two CQRs as well. 

The present paper is organized as follows. In Sect. 
II we present our theoretical model within the frame- 
work of CSDFT. In Sect. Ill we study the ground state 
properties of the CQR molecules. The phase diagrams of 
the ground state configurations of few electron quantum 
ring molecules in external magnetic field are obtained. In 
Sect. IV, we discuss the persistent current in the CQRs 
and we conclude our work in Sect. V. 



II. THEORETICAL MODEL 



our results qualitatively. The intra-ring and inter-ring 
Hartree potentials are given by 
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respectively, with the inter-ring distance d = |d|. 

In CSDFT all the quantities are functionals depending 
on the spin-up (p') and spin-down (p^-) densities, and the 

vorticity V(r) = £^-V x Therefore the exchange- 

correlation scalar and vector potential can also be writ- 
ten as a functional of these quantities, which are given 
respectively by 



Within the current spin-density functional theory 2 ^, 
we study the magnetic field dependence of the ground 
state of two vertically coupled GaAs quantum rings con- 
taining few electrons. The lateral electron confinement 
in the CQRs are described by the displaced parabolic 
potential model V(r) = \rn*uj^(r — ro) 2 in the xy-planc, 
where r = (x, y) = (r, 6), loq is the confinement frequency 
and ro is the radius of the ring. The two stacked identi- 
cal rings are coupled in the z direction with the potential 
V(z) described by two coupled symmetric quantum wells 
with a barrier of finite height. The quantum wells are 
assumed to be W = 120 A wide with the height of the 
barrier Vq = 250 meV between them. For these parame- 
ters we findi the following expression for the energy split- 
ting A = 22.86 exp[-d(A )/13.455] meV between the two 
lowest levels in the coupled quantum wells separated by 
a distance d. We also consider a homogeneous magnetic 
field B = Be z applied perpendicularly to the xy-planc, 
which is described by the vector potential A = Breg/2 
taken in the symmetric gauge. The Kohn-Sham orbitals 
ipjnm<j(j) = exp(— imd)(j) nmrT (r)Zj(z) are used to express 
the density and ground state energy. The Kohn-Sham 
equation in CSDFT for the CQRs is given by 
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where a =1 or J, is the z component of the electron spin 
and Lo c = eB/cm* is the cyclotron frequency. The total 
density in the rings is p(r) = Y, a J2n°m \4>nm<j{r)\ 2 . Be- 
cause we are adopting two identical rings, the density in 
each ring is half this total density. In the calculation, we 
approximate the density in the z direction by (5-functions 
at the center of the quantum wells. This approxima- 
tion has been used in our previous workl for two coupled 
quantum rings (QRs) with B — and will not change 
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where j p (r) is the paramagnetic current. Also in CS- 
DFT, an expression for the exchange-correlation en- 
ergy functional must be found and here we adopt the 
local density approximation remembering that in the 
bulk the total current density is zero (V = B). Thus 
the local density approximation can be implemented, 
which reads E xc = J dr^^-e xc [p^ , p^,V], where e xc is 
the exchange-correlation energy per particle of the uni- 
form two-dimensional electron gas in a magnetic field 
V = B. Following Refi 2 ^ we assume that e xc [p^ ,p^, \V\ = 
B] = {^Y M \pA + v^\p\p^\) /(l + ^ 4 ), where v = 
2irhcp/eB is the filling factor. The expression for e xc 



connects the fitted form of Levesque et alJ^ e 
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which is valid for large magnetic field, to the form given 
by Tanatar-Ceperley 2 ^ [p^ , p^\ valid for zero magnetic 
field. We expand the eigenfunctions ^> nmcr (r) in the well- 
known Fock-Darwin basis in order to solve the Kohn- 
Sham equation. 

The ground state (GS) energy of the coupled rings as 
a function of magnetic field is obtained from 

E(N) = J2 4™ -\J dvp(r) [Vk nt ™(r) + V^(r)] 
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and the paramagnetic current density is given by 
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FIG. 1: (Color online) Electron energy levels e„ |m of a single 
quantum ring as a function of the magnetic field. The ring 
radius is fixed at ro = 2oo- 

The measurable current density is equal to j(r) = j p (r) + 
(e/cm*)p(r)A(r), where the second part corresponds to 
the diamagnetic current density. 

In the z direction we consider only the two lowest lev- 
els of the quantum wells that connect the two quantum 
rings. They are the symmetric bonding level and the 
antisymmetric antibonding level. The contribution from 
excited states due to confinement in the z direction is 
neglected because the confinement in the z direction is 
much stronger than that in the plane. Therefore the mo- 
tion in the z direction may be assumed to be decoupled 
from the in-plane motion and the Kohn-Sham equations 
can be solved separately from the Schrodinger equation 
in the z direction that describes the tunneling between 
the quantum rings. In the limits of small and large inter- 
ring distance d, the results for single quantum rings arc 
recovered. On the other hand, in the limit of small ring 
radius (r — * 0), results of the CQDs are recovered^, 
too. 



III. GROUND STATE CONFIGURATIONS 

The energy levels e n ,m of a single quantum ring with 
fixed radius ro =2 ao (ro =300 A) are shown in Fig. 1 as 
a function of the magnetic field, where ao = yJh/m*ujQ. 
The energy levels are labeled by the radial quantum num- 
ber n — 0,1,2,..., and the angular quantum number 
m = 0, ±1,±2, .... The applied magnetic field breaks 
the ±77i degeneracy and leads to angular momentum 
transitions in the ground state. Although we assume 
a fixed value for the ring radius ro =2 ao, the effects 
observed in Fig. 1 remains the same for different ring ra- 
dius. The only difference is that the value of the magnetic 
field where the crossing occurs is rescalcd and assumes a 



smaller value when the ring radius is increased. Using 
the single-particle picture one expects already that the 
few electron ground state of the CQRs will be strongly 
affected by the magnetic field. In contrast, qualitative 
different behavior from two vertically coupled quantum 
dots was found only when iV > 8 in the absence of the 
magnetic field^. 

For the CQRs, the ground state configuration changes 
as a function of the inter-ring distance (tunneling en- 
ergy) and magnetic field. The phase diagrams presented 
in Fig. 2 show the different configurations of the ground 
state of CQRs for a fixed ring radius ro=2 ao and con- 
finement energy Hloq =5 meV as a function of the mag- 
netic field and the inter-ring distance for (a) N = 3, (b) 
N = 4, (c) N = 5 and (d) N = 6. To perform the nu- 
merical calculation, we use typical GaAs values for the 
effective mass m* = 0.067mo and the dielectric constant 
e s =12.4. For small magnetic field the GS configurations 
are the same as found previously for zero magnetic field^. 
With increasing magnetic field, many transitions in the 
ground state are observed, as can be seen in Figs. 2(a-d). 
The ground state phases are labeled by three quantum 
numbers (S Zl M z ,I z ): total spin S z , total angular mo- 
mentum M z and the isospin quantum number I z , which 
is the difference between the number of electrons in the 
bonding state and in the antibonding state divided by 2. 

For large inter-ring distance the two rings become de- 
coupled. On the other hand, when the distance between 
them is small, they are strongly coupled acting as a single 
one with isospin number I z = N/2 where all the electrons 
are situated in the bonding state. In the latter regime, 
six different ground state configurations with I z = 3/2 
arc found for 7V=3 (B < 4 T) in Fig. 2(a). Qualitatively 
the different GS configurations as a function of the mag- 
netic field can be understood through the single-particle 
(SP) picture in Fig. 1 together with Hund's rules. For 
example, when B < 1.37 T the SP energy level eo.o is 
lower than eo,i, hence it can be filled with two electrons 
with opposite spin and the third electron occupies the 
state eo,i with spin up, resulting in the GS configuration 
(1/2,1,3/2). But when B > 1.37 T, the energy e ,i is 
lower than eo.o an d it becomes energetically more favor- 
able to fill the state eo,i with two electrons instead of eo,o 
and the transition (1/2, 1, 3/2) -> (1/2, 2, 3/2) in the GS 
takes place. We also notice that in Fig. 2(a) this transi- 
tion occurs effectively at lower magnetic field B = 0.85 
T because of the many-body effects. The other tran- 
sitions (1/2,2,3/2) -> (1/2,4,3/2) -> (1/2,5,3/2) -> 
(1/2,7,3/2) -> (1/2,8,3/2) are due to the same mech- 
anism just explained. When the inter-ring distance in- 
creases, the difference between the bonding and anti- 
bonding states decreases and transitions in the GS con- 
figurations can be observed too. For large inter-ring 
distance in Fig. 2(a), the GS configurations with one 
electron in the antibonding state (I z = 1/2) are found. 
When the bonding-antibonding energy splitting A is less 
than the difference between two subsequent lateral bond- 
ing states the electron changes to the lowest unoccupied 
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FIG. 2: The phase diagram of the ground state of CQRs for a fixed ring radius (ro = 2ao) for: (a) N = 3, (b) N = 4, (c) N = 
5, and (d) N = 6. The insets indicate the three quantum numbers (S z , M z , I z ) which designate a ground state configuration. 



antibonding state yielding a GS transition. The state 
(3/2, 1, 1/2) can not be explained only using the SP pic- 
ture, because it has three aligned spins in different SP 
states. This is a clear manifestation of the many-body 
effects, where the total energy is minimized by the ex- 
change interaction, when two electrons are in the same 
quantum state. 

For N=A in Fig. 2(b) we observe eight different GS 
configurations for B < 4 T in the strong coupling regime 
(I z = 2) consistent with the SP picture and Hund's rule. 
The total spin in the z direction of these GS configu- 
rations alternates between and 1 as a function of the 
external magnetic field. Increasing the inter-ring distance 
we found GS configurations with one (J 2 = 1) and two 
(I z = 0) electrons occupying the antibonding sates. No- 
tice that for a fixed isospin I z =1 or 2, the total angular 
momentum changes by N/2 —2 as function of the mag- 
netic field. This fact is related to the even number of 
electrons (N =4) occupying the CQRs. In the strong 



coupling regime, the even number of electrons can be ar- 
ranged in such a way that S z oscillates between and 1 
and M z changes by N/2 with increasing the applied mag- 
netic field. Also the state (2, 2, 0) has all spins aligned 
due to the exchange interaction. 

For N=5 seven different GS configurations in the 
strong coupling regime (I z = 5/2) are found in Fig. 2(c) 
for B < 4 T. All of them are of total spin S z = 1/2. In 
this regime with increasing magnetic field, the total angu- 
lar momentum alternately changes by 3 [(1/2, 1,5/2) — > 
(1/2,4,5/2)] and by 2 [(1/2,4,5/2) (1/2,6,5/2)]. 
Since the number of electrons is =5, two SP-levels 
are filled by two pair electrons and the third level by a 
single electron, which causes a change of the angular mo- 
mentum by either 2 or 3 with increasing magnetic field. 
Increasing the inter-ring distance we found GS configu- 
rations with one (I z = 3/2) and two (I z = 1/2) electrons 
occupying the antibonding states. For these configura- 
tions the total spin can achieve values larger than 1/2. 
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FIG. 3: (Color online) M z /B as a function of magnetic field 
for different inter-ring separations d =30 A (red dash curves) , 
50 A (blue dotted curves), and 70 A (black solid curves). 



FIG. 5: The persistent current as a function of the magnetic 
field in the strong coupling regime (d = 30 A) for CQRs with 
different number of electrons: (a) N = 3, (b) N = 4, (c) 
N = 5, and (d) N = 6. 



The states with total spin S z = 3/2 and isospin I z = 3/2 
or I z = 1/2 are induced by the exchange interaction. 
Also the exchange interaction leads to the emergence of 
the maximum spin polarized state (5/2,4,1/2). For in- 
termediate (i-values, we have AM Z =1 transitions with in- 
creasing magnetic field corresponding to a single angular 
momentum increase of a single electron in this molecular 
state. 

For N = 6 in Fig. 2(d) nine different GS configura- 
tions arc found in the strong coupling regime (J z =3) 
for B < 4 T. Notice that for I z —3 the total angular 
momentum changes by N/2=3 as function of the mag- 
netic field due to the even number of electrons (iV=6) in 



CO 



CO 



N=3 



(a) 



N=5 



(c) 



N=4 



(b) 



N=6 



(d) 



12 3 

B (Tesla) 



1 2 3 

B (Tesla) 



FIG. 4: (Color online) The same as Fig. 3 but now for spin 

s.. 



the CQRs. Increasing the inter-ring distance we found 
the GS configurations with one (I z = 2), two (I z = 1) 
and three (I, = 0) electrons occupying the antibonding 
states. When I z = 2, the total spin is S z = 1 because the 
bonding states are occupied by 5 electrons and the GS 
configurations are the same as found earlier for JV = 5 in 
the strong coupling regime, which always have one elec- 
tron unpaired and the lowest antibonding state is filled 
with one more spin 1/2 electron. Many different states 
are induced by the exchange interaction when N — 6, 
e.g.,, the eight states composed with total spin 5*2 = 1 
or S z = 2 and isospin I z = 1. The GS configurations 
(3,6,0) and (3, 12,0) have all spins aligned. 

In phase (3,6,0), 3 electrons both in the bonding and 
the antibonding states occupy successive angular mo- 
mentum states and the total angular momentum M z = 
^(^ — 1) which is the densest spin polarized electron con- 
figuration available in a quantum ring molecule of 6 elec- 
trons. Such a state is referred to as the maximum den- 
sity droplet and was observed experimentally in a quan- 
tum dot in the presence of a magnetic fields However, 
we notice that there is no corresponding single quantum 
ring phase for such a configuration. It exists only in the 
CQRs because of the many-body interactions combined 
with the inter-ring quantum tunneling effect. In fact, the 
phases (5/2,4,1/2) for N = 5 found in Fig. 2(c), (2,2,0) 
for N = 4 in Fig. 2(b), and (3/2,1,1/2) for N = 3 in 
Fig. 2(a) are the maximum density droplet states in the 
QR molecules. 

In order to show more clearly the behavior of the total 
angular momentum and total spin in the CQRs, we plot 
M z /B and S z as a function of magnetic field in Fig. 3 
and Fig. 4, respectively, for different inter-ring distances 
d =30 A (the dash curves), 50 A (the dotted curves), 
and 70 A (the solid curves). The oscillation of M z /B 
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FIG. 6: (Color online) The persistent current (solid curves) as 
a function of the magnetic field in the weak coupling regime 
(d = 70 A) for CQRs with different number of electrons: (a) 
N = 3, (b) N = 4, (c) N = 5, and (d) N = 6. The dashed 
(red) and dotted (blue) curves correspond to the bonding and 
anti-bonding currents, respectively. 



in Fig. 3 is related to the magnetization of the CQRs 
and the corresponding persistent current which will be 
shown below. In Fig. 4 we see that with increasing the 
distance between the two rings the spin polarization of 
the system is enhanced. Generally, the QR molecule in 
the weak coupling regime (d =70 A) is spin-polarized in 
a wider range of magnetic field than that in the strong 
coupling regime (d =30 A, which is practically in the 
atomic regime). The maximum spin polarization S z = 
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FIG. 7: (Color online) The persistent current obtained 
through CSDFT calculation (solid black curve) and the single 
particle approximation (dashed green curve), considering the 
CQRs with four electrons in both (a) strong (d = 30 A) and 
(b) weak coupling regime (d — 70 A). 
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FIG. 8: The total spin-current considering (a) N = 3 in the 
strong coupling regime (d = 30 A) and (b) for TV = 6 in the 
weak coupling regime (d = 70 A). 

N/2 always occurs in the molecular phases in the weak 
coupling regime at finite magnetic field. 



IV. PERSISTENT CURRENT 

In Fig. 5 the persistent current is shown, determined 
by^S I = J j{ r )dr, as a function of the magnetic field 
in the strong coupling regime (d = 30 A) for different 
number of electrons N = 3, 4, 5 and 6. The diamagnetic 
contribution to the persistent current in our model can be 
evaluated analytically and is given by I d = £^- c . When 
the number of electrons in the CQRs is odd (N=3, 5), the 
persistent current oscillates with linear segments, and the 
different segments appear because of the change in the 
total angular momentum. As can be viewed in Figs. 3(a) 
and 3(c), for d = 30 A the total angular momentum 
changes as function of the magnetic field. The total mo- 
mentum increases from Mf to M> with increasing mag- 
netic field following the rule M> = M< + (N ± l)/2 
indicated in Fig. 2. These relations for the total angu- 
lar momentum in the strong coupling regime can be ex- 
plained by the breaking of the ±m degeneracy and the 
crossing between states with angular quantum number 
(m + 1) and ±m. For example, in the case of N = 5 and 
B < 0.8 T the SP energy levels can be filled in the follow- 
ing way: two electrons in the state eo,0; two in eg,i and 
one in eo.-i, which gives the GS (1/2, 1, 5/2). Increasing 
the magnetic field the state eo,2 crosses the state eo,-i 
at B = 0.8 T and now the configuration of the GS is 
(1/2,4,5/2) with: two electrons in the state erj,0j two in 
eo,i and one in eo,2- For B > 1.37 T the level eo,2 crosses 
the eo,i and the configurations of the GS is (1/2,6,5/2) 
with: one electron in the state eo,Oj two in eo.i and two in 
e.Q.2- Therefore the total angular momentum jumps from 
M z = 1 to M z = 4 and from M z = 4 to M z = 6 respect- 
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FIG. 9: (Color online) The persistent current for CQRs with: 
(a) TV = 3, (b) N = 4, (c) TV = 5, and (d) TV = 6, as a 
function of the distance between the two coupled quantum 
rings. Also, we consider two different values of magnetic field 
B — IT (black solid curve) and B = 1.1 T (B = 1.25 T) (red 
dashed curve) in the upper (bottom) panel. 



ing the crossing between the SP energy levels. However, 
the value of the magnetic field at which the GS transi- 
tion occurs is not exactly the same as the SP crossings 
because it is reduced by the electron-electron interaction. 

When the number of electrons is even TV = 4 [see 
Fig. 5(b)] or TV = 6 [see Fig. 5(d)] in the strong coupling 




20 30 40 50 60 70 30 40 50 60 70 



d(A) 



d(A) 



FIG. 10: (Color online) The persistent spin-current for CQRs 
with: (a) TV = 3, (b) TV = 4, (c) TV = 5, and (d) TV = 6, as 
a function of the distance between the two coupled quantum 
rings. Also, we consider two different values of magnetic field 
B = 1T (black solid curve) and B = 1.1 T (73 = 1.25 T) (red 
dashed curve) in the upper (bottom) panel. 



regime (d = 30 A), the persistent current also oscillates 
with linear segments, but an additional fine structure is 
found. In this case, the difference between two consecu- 
tive total angular momentum transitions as a function of 
the magnetic field is always equal to M> — Mf = TV/2 
as shown in Figs. 5(b) and 5(d). When the total spin 
is different from zero (S z ^ 0), the ground state is un- 
stable because it corresponds to an open-shell system, 
and consequently an intermediate structure in the total 
persistent current appears. For TV-even, S z oscillates be- 
tween and 1 [see the dash curves in Figs. 4(b) and 4(d)] 
and these transitions are reflected in the total persistent 
current. The more stable linear segments correspond to 
closed-shell configurations (S z = 0). 

The persistent current as a function of the magnetic 
field in the weak coupling regime (d = 70 A) for differ- 
ent number of electrons TV = 3, 4, 5 and 6 are shown 
in Figs. 6(a-d), respectively. The dotted and dashed 
curves correspond to the bonding and antibonding cur- 
rents, respectively. The bonding (antibonding) current 
is the contribution of the electrons in the bonding (an- 
tibonding) state to the persistent current. For example, 
the bonding paramagnetic current density is found by 
considering only the contribution of the bonding states 
in Eq. ([5]), and so on. The solid curve represents the total 
current, which is the sum of the bonding and antibond- 
ing currents. In the weak coupling regime, the bonding 
and antibonding states can be occupied simultaneously 
as can be viewed through the changes in I z indicated in 
Figs. 2(a-d). Therefore, different combinations of the to- 
tal angular momentum are possible as a function of the 
magnetic field and any change of M z leads to a jump 
in the persistent current. For large values of the mag- 
netic field, such jumps become more frequent because 
of a rapid variation of the angular momentum at large 
magnetic fields and consequently the amplitude of the 
oscillation in the the persistent current is reduced. 

In order to understand the relevance of the many- 
body effects, wc show in Fig. 7 the persistent current 
calculated using CSDFT (solid curve) and using the sin- 
gle particle approximation (dashed curve), considering 
the CQRs with four electrons in both the strong [see 
Fig. 7(a)] and the weak [see Fig. 7(b)] coupling regime. 
When the single-particle approximation is considered, we 
note that most of the little jumps in the persistent cur- 
rent are missing and in some regions the single-particle 
results give an opposite persistent current compared to 
the persistent current evaluated by CSDFT. Therefore, 
through this comparison in Fig. 7 wc clearly sec the im- 
portance of the many-body interactions when estimating 
the persistent currents in CQRs. 

We present in Fig. 8 the persistent spin-current 
for a CQRs with (a) TV = 3 in the strong coupling 
regime (d = 30 A) and (b) TV = 6 in the weak cou- 
pling regime (d = 70 A). The spin-current is defined 
as the difference between the current of electrons 
with spin-up and with spin-down. The paramagnetic 
spin-current density is given explicitly by j^(r) = 
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The diamagnctic contribution for the persistent spin- 
current is given by /J = eB ^ m ~^ V . The persistent 
spin-current is zero in the magnetic field ranges where 
the total spin S z in a certain GS configuration is zero 
as can be viewed in Fig. 8(b). When the CQRs is 
occupied by an odd number of electrons, the persistent 
spin-current is never zero [see Fig. 8(a)] because now 
always a not-aligned electron is present. 

In addition to the external magnetic field, the dis- 
tance between the two quantum rings is another impor- 
tant parameter affecting the persistent current and the 
spin-current. Fig. 9(a) shows the persistent current for 
a CQRs with N = 3 for different values of the mag- 
netic field (a) B = 1 T (solid curve) and B = 1.1 T 
(dashed curve). The respective spin-currents are pre- 
sented in Fig. 10(a). For the CQRs with N = 4 we 
plot the same quantities in Fig. 9(b) and Fig. 10(b) for 
the magnetic fields B = 1 T (solid curves) and B = 1.1 
T (dashed curves). We plot the persistent current (spin- 
current) in Fig. 9(c) (Fig. 10(c)) for a CQRs with N = 5 
for the magnetic fields B=l T (the solid curves) and for 
B = 1.25 T (the dashed curves). The persistent current 
in Fig. 9(d) and persistent spin-current in Fig. 10(d) are 
evaluated for different values of the magnetic field B=l 
T (solid curves) and B = 1.25 T (dashed curves) for 
a CQRs with N = 6. In Figs. 9-10 we notice that for 
a fixed value of the magnetic field, the persistent cur- 
rent (spin-current) as a function of the inter-ring dis- 
tance is practically constant and exhibits a jump when 
a ground state transition occurs. Therefore, by vary- 
ing the inter-ring distance for fixed value of the applied 
magnetic field, the GS configuration of the CQRs can 
be determined. For N = 3 and B = 1 T (B = 1.1 
T), we have the transition (1/2,2,3/2) -> (3/2,1,1/2) 
((1/2,2,3/2) -> (1/2,3,1/2)) as a function of the inter- 
ring distance. When B = 1 T the total spin changes 
from S z — 1/2 to S z — 3/2 and the total angular mo- 
mentum from M z = 2 to M z = 1. Therefore the total 
angular momentum decreases and the persistent current 
increases [see Fig. 9(a)], because of the reduction in the 
paramagnetic current [see Eq. ©]. The increase in the 
total spin in the z-direction causes an amplification of the 
persistent spin-current, as can be viewed in Fig. 10(a). 
On the other hand, when B = 1.1 T the total spin does 
not change and the total angular momentum decreases, 
which leads to a reduction of both the persistent current 
and the spin-current. 

When the CQRs are filled with 4 electrons and the ap- 
plied magnetic field is B — 1 T, the transition (0, 2, 2) — > 
(2,2,0) occurs as a function of the inter-ring distance. 
Increasing the magnetic field slightly (B = 1.1 T), an in- 
termediate GS configuration appears and now we have 
the transitions (0,2,2) -» (1,3,1) -> (0,4,0). When 
B = 1 T, the total spin changes from S z — to 
S z = 2 and the total angular momentum is not af- 
fected. Thus the persistent current is almost unaltered 



[see Fig. 9(b)] and the persistent spin-current changes 
abruptly [see Fig. 10(b)], because all the carriers are 
spin polarized for S z = 2. Moreover, when B = 1.1 T 
the persistent current firstly decreases (56A< d <65A) 
and afterwards increases id >65A) due to the increase 
of the angular momentum. The persistent spin-current 
is zero for d <56A, because the net spin is zero in this 
range. For 56A< d <65A the persistent spin-current 
is different from zero and increases more for d >65A, 
because of full alignment of the electron spins. For 
N = 5, we choose two different magnetic fields B = 1 
T and B = 1.25 T which gives the following GS tran- 
sitions (1/2,4,5/2) -> (1/2,3,3/2) -> (5/2,4,1/2) and 
(1/2,4,5/2) -> (3/2,5,3/2) -> (5/2,4,1/2), respectively. 
The initial and final GS configurations are the same for 
these two magnetic fields and both have the same to- 
tal angular momentum M z = 4, and only the GS con- 
figurations with I z = 3/2 are different. This difference 
brings an opposite behavior of the persistent current [see 
Fig. 9(c)] for the I z = 3/2 GS configurations and in one 
case (B=l T) it increases and the other (B = 1.25 T) de- 
creases, due to the difference of total angular momentum 
in both cases. The latter-type behavior of the persis- 
tent spin-current [see Fig. 10(c)] in both cases (B = 1 
T and B = 1.25 T) is a consequence of the increase 
in the spin polarization as a function of the inter-ring 
distance. The ground state transitions for the applied 
magnetic field B = 1 T and B = 1.25 T are respec- 
tively (0,6,3) -> (1,5,2) -> (2,5,1) -> (3,6,0) and 
(0, 6, 3) -> (1, 7, 2) -> (2, 7, 1) -> (3, 6, 0), when the CQRs 
are occupied with N = 6 electrons. Again the initial and 
final GS configurations for the magnetic fields chosen are 
the same and in this case two intermediate configura- 
tions appear. When B = 1 T (B = 1.25 T), the persis- 
tent current [see Fig. 9(d)] increases (decreases) when the 
inter-ring distance corresponds to one of the intermediate 
configurations, because of the increase (decrease) in an- 
gular momentum. The persistent spin-current is zero for 
the GS configuration (0, 6, 3) due to the spin-unpolarizcd 
electrons [see Fig. 10(d)]. The persistent spin-current 
[see Fig. 10(d)] of the intermediate GS configurations 
(1,5,2) and (2,7,1) ((1,7,2) and (2,5,1)) decreases (in- 
creases) because the paramagnetic spin-current is larger 
(smaller) than the diamagnctic spin-current. The GS- 
configuration (3,6,0) for B = 1 T has a small spin- 
current due to the compensation of the paramagnetic 
spin-current by the diamagnctic spin-current, what does 
not happen for B = 1.25 T, where the diamagnctic spin- 
current is larger than the paramagnetic spin-current. 



V. SUMMARY 

We studied systematically the electronic structure and 
the persistent current in two vertically coupled quantum 
rings as a function of an applied external magnetic field 
and the inter-ring distance. A rich variety of ground state 
configurations were found, that are induced by changing 
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the inter-ring distance and the applied magnetic field. 
We found that these transitions are governed by some 
rules, which are related to the even/odd number of elec- 
trons, the single-particle picture, Hund's rules and many- 
body effects. The found ground state configurations are 
summarized in phase diagrams, which generalize the pre- 
vious B=0 results^ to non-zero values of the magnetic 
field. Moreover, the persistent current and spin-current 
were evaluated for the CQRs and we showed that their 
variation with magnetic field is governed by the values of 
the total angular momentum and total spin. Therefore, 
once the GS configuration is known the persistent cur- 
rent and spin-current can be estimated. Furthermore, 
we provide useful information to determine the ground 
state configuration through measurements of the persis- 
tent current. Also we showed that the electron-electron 
interaction strongly influences the sign and size of the 
persistent current. 

The simple rules and the complete set of results for 



iV=3,4,5, and 6 electrons shown in the present work may 
not be compared to available experimental data. The rea- 
son for that is because we assumed a system composed 
of two symmetric QRs vertically coupled containing few 
electrons and until now the experimental fabrication of 
such an arrangement is quite difficult. Nonetheless, we 
believe that our results and conclusions may be very use- 
ful to interpret experimental data when those difficulties 
are overcome. 
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